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New theoretical results from large-scale relativistic close coupling calculations reveal the pre- 
cise effect of resonances in collisional excitation of x-ray lines of Ne-like Fe xvn. Employing the 
Breit-Pauli R-matrix method and a 89-level eigenfunction expansion including up to n = 4 levels 
shows significant resonance enhancement of the collision strengths of forbidden and intercombina- 
tion transitions. The present results differ from all previous calculations, heretofore without detailed 
resonance structures, and should help resolve longstanding discrepancies. In particular, the present 
line ratios of three benchmark diagnostic lines 3C, 3D, and 3E at 15.014, 15.265, and 15.456 A 
respectively, are in excellent agreement with two independent measurements on Electron-Beam- 
Ion- Traps [Laming etal. , Astrophys. J 545, L161 (2000) and Brown etal. , Astrophys. J 502, 
1015 (1998)]. The strong energy dependence due to resonances in these and other cross sections is 
demonstrated for the first time. It is of general importance and strongly manifests itself in x-ray 
plasma diagnostics. 

PACS number(s): 34.80.Kw 

Prominent Fe xvn X-ray lines have long been observed from laboratory and astrophysical sources [e.g. [[|-[|], most 
recently in the Chandra and the XMM-Newton spectra of active galactic nuclei, x-ray binaries, supernovae, and flar- 
ing and non-flaring active stellar coronae jsjQ] - In spite of a number of theoretical, experimental, and observational 
studies over decades there remain outstanding discrepancies related to atomic processes and astrophysical effects 
responsible for Fe xvn spectral formation in high temperature sources [1-15]. However, the interpretation of experi- 
mental and observational ratios of line intensities usually relies on collisional-radiative (C-R) models using theoretical 
cross sections that neglect the fundamental role of resonant excitation, which preferentially affects the forbidden and 
intercombination transitions as opposed to dipole allowed ones. Of particular interest are three prominent x-ray tran- 
sitions to the ground level ls 2 2s 2 2p 6 % from excited levels: 3C (A 15.0141) ls 2 2s 2 2p 5 [l/2]3d 3/2 X P? (level 27), 3D 
A 15.2651 ls 2 2s 2 2p 5 [3/2]3d 5/2 3 D? (level 23), and 3E A 15.4561: ls 2 2s 2 2p 5 [3/2]3d 5/2 3 P? (level 17), (level numbers 
denote energy order relative to the ground level 1). While the 3C is dipole allowed, the 3D and 3E are spin-forbidden 
intercombination transitions. All but one of the previous calculations used the distorted wave (DW) approximation 
that neglects channel coupling, and hence resonances Jq— |Tl[] ■ A point of confusion has been that the only previous 
coupled channel (CC) calculation [jl2| also yielded the same results as the DW calculations |l]{] , because the CC 
R-matrix calculations were at energies above all target levels in the eigenfunction expansion; and therefore resonances 
were not included. Since both the DW and the CC results agreed, it has hitherto been assumed that resonance effects 
are not important, and that the cross sections are relatively constant with energy (e.g. P,|l3|]). On the other hand, the 
experiments on Electron-Beam- Ion- Traps (EBIT) at Lawrence Livermore National Laboratory (e.g. Q), and at the 
National Institute of Standards and Technology (NIST, 0), do not measure absolute excitation cross sections directly. 
Rather, they measure relative intensity ratios R1=3C/3D and R2=3E/3C at a few selected energies. Therefore the 
presence of resonances, and more generally the energy variations, are not readily discernible in experimental data 
(discussed later). The measured values differ from theoretical ones by up to 50% for Rl, and a factor of 2 for R2 
|p],p| Jl^JT4] ,|6|| . In order to accout for these discrepancies several mechanisms based on atomic and astrophysical effects 
have been put forward, such as polarization, resonance scattering, and dielectronic satellite blending, which may be 
of varying importance under appropriate conditions [p]J^j7|,|l5]] . 

However, as evident from the results in this Letter, there are extensive resonance structures in cross sections 
for all transitions in Fe xvn due to many Rydberg series of resonances converging on to a number of n = 3 and 
ii = 4 levels. Many infinite and interacting series of resonances arise due to coupling between open and closed 
scattering channels which, in principle, must be included in order to obtain the cross sections precisely. But for 
a complex ion such as Fe XVII the number of channels is very large and the CC calculations become enormously 
difficult, especially since relativistic fine structure must also be considered in addition to other atomic effects. Using 
the Breit-Pauli R-matrix (BPRM) method fll|[ijj we construct a large eigenfunction expansion including 89 levels 
corresponding to 49 LS terms up to the n = 3 and the n = 4 complexes of Fe XVII. Full details of the calculations 
will be presented elsewhere, but we briefly outline these below. The coupled-channel wavefunction expansion for the 
(e + Fe XVII) system may be expressed as $(E; e + Fe xvn) = £\ Xi(Fe xvn)# c (£j) + Ylj c j^j( Fe XVI )> where the 
'F denote the continuum (E > 0) states of given each total angular momemtum and parity lit, expanded in terms 
of the core ion eigenfunctions Xii^iLiJi)', the 9 e (£i) refer to the free-electron partial wave, and the <E>j are short- 
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range correlation functions that also serve to compensate for orthogonality constraints. The 89 levels belong to the 
configurations 2s 2 2p 6 , 2s 2 2p 5 (3s, 3p, 3d), 2s 2 2p 5 (4s, Ap, Ad, 4/), 2s 1 2p 6 (3s, 3p, 3d), 2s 1 2p 6 (4s, 4p, 4d, 4/). We consider 
total symmetries 2J< 51 of both parities explicitly in the BPRM calculations. The size of the calculations (possibly 
the largest (e + ion) scattering calculations to date) may be gauged by the fact that the dimension of the Hamiltonian 
matrices ranges up to 10286, for J = 3.5 with 395 free channels and 486 bound channels (<F,) (the largest number of 
free channels per symmetry is 401); 25 continuum basis functions are used to represent the ^(e + Fe XVll) in the 
inner R-matrix region. Relativistic distorted wave (RDW) and Coulomb-Born-Bcthe approximations are employed to 
'top-up' the partial wave summations. Particular attention is paid to the resolution of resonances, with cross sections 
computed at about 20,000 energies. 

Figs, l(a-c) display the dense resonance structures in excitation collision strengths for the 3C, 3D, and 3E transitions 
respectively. The results are compared with the RDW values (filled circles, that essentially represent previous DW 
calculations (such as from the HULLAC code ]9],[l0|), and the one non-resonant value from the previous CC R-matrix 
calculation (filled square, [0). As seen from the figures, the other data correspond roughly to the 'background' 
collision strength compared to the present detailed results. The energy behavior of the three transitions may be 
ascertained from the eigenfunction expansions of the upper levels: 3C: 0.7857|27) + 0.1753|23) + 0.0305|17), 3D: 
0.7479|23) + 0.2010|27) + 0.0491|17), and 3E: 0.9150|17) + 0.0767|23) + 0.0030|27), where the eigenkets correspond to 
the energy level indices. As the mixing coefficients indicate the 3D has a significant component of ^P ), but the 3E 
considerably less so. Therefore the 3D would depart from LS coupling to intermediate coupling to jj-coupling schemes 
along the neon isoelectronic sequence with Z, and with f2(3D) increasing with energy like f2(3C). On the other hand 
the 3E remains largely a (spin) forbidden transition and Q(3E) decreases with energy. The Einstein A-coefficients of 
the 3C, 3D, and 3E are 2.47 x 10 13 , 6.01 x 10 12 , and 9.42 x 10 10 sec -1 respectively. Other interesting atomic physics 
aspects of these and other transitions in neon-like ions will be discussed in later publications, with particular reference 
to laser transitions. 

Spectral formation in laboratory and astrophysical plasmas needs to be distinguished from each other. While we 
describe the latter using a C-R model for Fe XVII, generally using collision strengths averaged over a Maxwellian 
characterizing a temperature, the former are measured using a mono-energetic beam with a certain veocity spread 
assumed to be a gaussian. We compute the collision strengths averaged three ways: (i) Maxwellian, (ii) gaussian, and 
(iii) numerical. Present results for (ii) and (iii) are shown in Figs, l(a-c). In order to compare with EBIT experiments 
the gaussian FWHM is taken to be 30 eV. We note that the energy variations, dependent on resolution and density 
of resonances in a given energy region, are reflected in the gaussian averages (ii); the numerical averages are slowly 
varying. The theoretical line ratios using gaussian and numerical averages (a ■ v) (GA and NA respectively; a denotes 
the cross section), are computed taking account of the radiative branching ratios for 3C, 3D, and 3E which are: 1.0, 1.0, 
and 0.89 respectively. Therefore we have Rl = 3C/3D = (a 3C ■ v)/(a 3D ■ v) and R2 = 3E/3C = 0.89(cr 3E ■ v)/(a 3C ■ v). 
The results are compared with several EBIT measurements in Table I. These direct theoretical results agree with 
measured values to within experimental errors. Although our calculations extend up to the n = 4 levels, the highest 
threshold (level 89) is at ^84.5 Ryd or 1.15 keV; no resonances above this energy are therefore included. While 
the line ratios at 0.85 keV and 0.9 keV are computed directly from the collision strengths as in Figs, l(a-c), some 
extrapolation of the averaged values is made to compare with the experimental value at 1.15 keV since we expect 
similar resonance enhancement above the n = 4 due to still higher thresholds. The GA results differ rather more 
than the NA results but still agree with experiment, in contrast to all other theoretical values. However, there is 
significant variation with energy in the GA collision strengths as they 'oscillate' irregularly with energy depending 
on the density of resonances prevalent in a given range. This appears to be reflected in similar 'oscillatory' structure 
seen in the experimental values reported in |T^| throughout the energy range 0.1 - 4 keV. 



TABLE I. Comparison of present line ratios for R1=3C/3D and R2=3E/3C with EBIT measurements 
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f E3|; ^present values with extrapolation of resonance enhancement from ab initio collision strengths from E < 1.02 
keV (see text). 

We also find that it is the n = 4 resonances, rather than the n = 3, that dominate the enhancement of electron 
impact cross sections in Fe XVII. The n — 4 resonances begin with resonant configurations 2s 2 2p 5 3£M' that manifest 
themselves from ^47 Ryd, considerably below the excitation thresholds of the 2p — > 3d lines 3C, 3D, and 3E at 
^60 Ryd. Therefore there arc relatively fewer n — 3 resonances and the n = 4 resonances greatly influence the 
near-threshold behavior of these cross sections. In the region 75 - 84.5 Ryd the resonances are due to thresholds 
corresponding to two-electron-excitation configurations 2s2p e A£ with weakly coupled channels and resonances. This 
region is expected to be dominated by n > 4 resonances. It might also be mentioned that radiation damping of 
autoionizing resonances in cross sections of Fe xvn is known to be negligible ]l8|-|2"cfl . 

We note a few specific cases that exemplify the findings in this Letter. An application of the present rates would be 
to estimate more preciesly the degree of resonance scattering of the 15.014 A line in the solar and stellar coronae as 
discussed in jj],[l3| . Based on our C-R model with the present collisional data and extensive new radiative calculations 
for Fe xvii transition probabilities using the BPRM method (as in and from the SUPERSTRUCTURE code 
p2fl , we obtain the 3C/3D ratio to vary between 2.63 - 3.20, which agrees with the values from EBIT and from flaring 
solar corona, but is still higher than the observed value jl5) of 1.87±0.21 from non-flaring active region. This further 
suggests that resonance scattering of 3C, or possibly some other mechanism, may be operative under astrophysical 
conditions . Another example is the discussion of the diagnostic utility of the 3C /3D line ratio by Brown et al. 
p3| who (a) assume the cross sections to be relatively constant, and (b) mutiply the RDW cross section of Zhang 
and Sampson by 1.25 to agree with the measured value. From the present work neither (a) nor (b) are necessary. 
Furthermore, Brown etal. obtain the temperature dependence of 3C/3D including an inner-shell satellite line blended 
with the 3D, but using constant cross sections. A revised analysis with the present cross sections should yield a 
different temperature diagnostics. 

More generally, the excitation of most of the n = 3 levels of Fe xvii is similarly affected. Using our 89-level 
C-R model, with newly computed collsional and radiative data employing the BPRM method, we investigated the 
prominent transitions corresponding to the 'coronal' x-ray lines 2s 2 2p 5 3s — ► 2s 2 2p e at A A 16.780, 17.055 and 17.100 A 
(e.g. f23]| ). Fig. 2(a) shows the detailed collision strength for the forbidden J = 1 — ► A 16.780 A (3F) transition, 
compared with previous DW values (filled circles Q and square |ll]]). The forbidden transitions are expected to 
be most enhanced by resonances, as in Fig. 2(a), since the background cross sections are much smaller than for 
allowed transitions. The temperature dependence of the forbidden (3F) to the allowed (3C) line ratio R3 = 3F/3C 
= 1(16. 780)/I(15. 014) is demonstrated in Fig. 2(b), and compared with that calculated using DW cross sections 
(filled squares pi|); at higher temperatures the resonance contribution decreases progressively with energy due to 
the Maxwellian factor. The electron density dependence is small; solid-line and dot-line correspond to 10 13 and 10 9 
cm~ 3 respectively. The 4 open circles with error bars are observed and experimental values. At all temperatures T 
< 10 K the present line ratio departs considerably from those using DW data without resonances, to more than a 
factor of 3 at about 10 6 K — a fact of considerable importance in photoionized X-ray plasmas that have temperatures 
of maximum abundance much lower than that in coronal equilibrium T m ~4x 10 6 K for Fe xvii, as marked. 



TABLE II. Rate coefficients at T c =200 eV for resonant excitation to the 2s 2 2p 5 3s levels 
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Resonant excitation may be considered indirectly by approximate methods such as quantum defect theory |E0[ or 
isolated resonance approximations (IRA) |4|]9|]ic}l . Table II compares the present rate coefficients for the individual 
3F line, and for all transitons to the 2s 2 2p 5 3s levels, with IRA calculations using HULLAC (|], multi-configuration 
Dirac-Fock jlo|], and scmi-relativistic Hartree-Fock Q| methods. The present data are ^50% higher and differ from 
P,fio[ which (i) neglect interference between resonances, and (ii) take limited account of the several infinite Rydberg 
series of resonances resulting in a considerable underestimate of resoance enhancement. The imprecision of the IRA 
is also indicated by more than a factor of 3 discrepancy between and j| , both using IRA but with different decay 
channels. Present results are about a factor of 2 lower than Although resonances are generally less important for 
highly charged ions, the resulting enhancement in the BP R-matrix cross sections for forbidden and intercombination 
lines, compared to the allowed line, yields the ratios Rl and R2 closer to experiments compared to other works. 
We also note that ionization and recombination processes are not included in the present C-R model; dielectronic 
recombination from Fe XVIII to Fe XVII may also contribute to line emissions. 

The main conclusions of this Letter are: (I) resonances are present in Fe xvn collision cross sections at all energies, 
(II) the energy dependence of different types of forbidden, intercombination, and allowed spectral transitions is 
demonstrated, (III) the effect on diagnostic line ratios is significant, and may be up to several factors compared 
to previous data neglecting resonance effects, (IV) the present work is expected to be generally applicable to the 
temeprature and density analysis of many astrophysical and laboratory X-ray sources (especially photoionized plasmas 
and forbidden lines) not only for Fe XVII but also other neon-sequence ions of importance in x-ray lasers |E 3 |, (V) 
this work is part of a large-scale program to compute excitation, photoionization, (c + ion) recombination |19| , po| , 
and transition probabilities of Fe xvn using the accurate BPRM method; for example, a set of A-values for ^30 000 
transitions have been computed. 
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FIG. 1. Breit-Pauli R-matrix collision strength f2 for 3C, 3D and 3E lines with detailed resonance structures as a function of 
incident electron energy; filled dots are RDW results. The gaussian (FWHM = 30 eV) and numerical averages are also shown 
(see text). 
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FIG. 2. (a): Breit-Pauli R-matrix collision strength Q for the 3F line; the filled circles and square are non-resonant DW 
calculations; (b): line ratio 3F/3C vs. T from a 89-level C-R model. The electron densities for solid-line and dot-line curves 
are 10 13 and 10 9 cm -3 respectively. The 4 open circles with error bars are observed and experimental values: from the solar 
corona at T m ~ 4MK ||, from the corona of solar-type star Capella at ~ 5MK |J, and from the EBIT experiment at 0.9 keV 
(log T = 7) 0. The filled squares are values using DW results 
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